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ABSTRACT 

Lithium is made up in the envelopes of massive Asymptotic Giant Branch (AGB) 
stars through the process of Hot Bottom Burning. In Globular Clusters, this process- 
ing is one possible source of the hot-CNO burning whose nuclear products are then 
ejected into the intracluster medium and take part in the formation of a second stel- 
lar generation, explaining the peculiar distribution of chemical elements among the 
cluster stars. We discuss the lithium yields from AGB stars in the mass range 3- 
~6.3M0, and from super-AGB stars of masses in the range 6.5-9 for metallicity 
Z=10~^. The qualitative behaviour of these yields is discussed in terms of the phys- 
ical structure of the different masses. Although many uncertainties affect the other 
yields of these stars (e.g. O, Na and Mg), even larger uncertainties affect the lithium 
yield, as it depends dramatically on the adopted description of mass loss. When we 
adopt our standard mass loss formulation, very large yields are obtained especially for 
the super-AGB stars, and we discuss their possible role on the lithium abundance of 
second generation stars in globular clusters. 

Key words: Stars: AGB and post-AGB; Stars: abundances; Globular clusters: gen- 
eral 



1 INTRODUCTION 

In highly evolved giants, the pristine lithium has been al- 
ready well depleted from the stellar atmosphere due to 
convective dilution with the internal layers where it has 
been burned. If lithium is observed, its abundance is as- 
cribed to temporary prod uction due to the action of the 
ICameron fc Fowled ([iQTJ) mechanism: ^Be is produced by 
fusion of ^He with *IIe, and rapidly transported to stellar re- 
gions where it can be converted into ^Li by electron-capture. 
Envelope models of asymptotic giant branch (AGB) stars 
ijScalo et al.lll975l ). in which the temperature at the bottom 
of their hot bottom convective envelopes (Tbco) reaches the 
hydrogen burning layers, show that the ^IIe(a,7)^Be can 
act (at Tbco~40MK). These models were able to explain 
the high lithium abundances found in some luminous red gi- 
ants, and the process took the name of Hot Bottom Burning 
(HBB). A differe nt, slow mixing pro cess named "cool bot- 
tom processing" (jNollett et al.|[2003l ) is instead at the basis 
of the hthium abundances seen in lower luminosity red gi- 
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ants (e.g. IWasserburg et all Il995l : ISackmann fc Boothrovdl 
Il999l) . The physical reasons of this process are not well 
studied, while the nucleosynthesis in HBB is simply based 
on time-dependent mixing in plain convection zones. Apart 
from being a key process for the lithium production, HBB is 
the key mechanism necessary to process to nitrogen the car- 
bon present in the envelopes of carbon stars, that becomes 
active at slightly larger Tbce~65MK. Observationally, lumi- 
nous AGB stars lose their carbon star character when they 
are lithium rich, although a small luminosity interval may be 
characterized by star s that are both lithium and carbon rich 
I Ventura et al.lll999h . as in fact occurs in the most luminous 
fraction of the J stars ijAbia et al.lll99ll : ISmith et al.lll995l ). 
In the most massive AGB stars, and especially in low metal- 
licity AGB stars, Tbcc becomes larger than ~80MK, and also 
the ON chain of the CNO cycle becomes active. In these 
envelopes, oxygen is cycled to nitrogen, and its abundance 
can be dramatically reduced. Thus some models attribute 
the existence of anomalous stars in globular clusters (GCs), 
having low oxygen and high sodium, to the formation of a 
second stellar ge neration (SG) in cluding matter processed 
by HBB (e.g. iVentura et al.|[200ll) . Other models attribute 
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Figure 1. Lithium abundances as a function of [Fc/H] in halo 
stars and in scarcely evolved stars in three GCs. Halo data are 
from Melenzez et al. (2009), represented as black triangles (LTE 
models) or grey (cyan) triangles (non LTE models). (Blue) full 
squares are from Sbordone et al. (2009), analyzed by 3D non LTE 
models. The two dashed horizontal lines represent an eye fit of 
the Melendez LTE (upper line) and non LTE model data (lower 
line) in the range of the GC abundances. Data for NGC 6397 arc 
from Lind et al. (2009). The three open triangles are values for 
three subgiants, and may not represent the turnoff abundances in 
this cluster. Data for NGC 6752 are from Pasquini et al. (2005), 
plotted as open or full squares according to the two different tem- 
perature scales used in this work. The full circles are the data for 
47 Tuc by Bonifacio et al. (2007). The limits of the lithium range 
in the 50 stars recently examined by D'Orazi et al. (2009) arc 
also given. Open circles represent the average abundances in the 
ejecta of models of 4, 5 an 6MQ(see Table 1) for three different 
mass loss rate formulations. 



the formation of the SG to the ejecta of fast ro tatinR mas- 
sive stars (FRMS, see e.g. lPecressin et al.ll2007ah . or even to 
pollution from gas expelled during hishly non-co nservative 
evolution of massive binaries (Ide Mink et al.l2009l ). although 
this latter model in particular can not explain the very high 
fracti on of SG stars present in most of the GCs so far exam- 
ined (jCarretta et al.ll2009al lbh. The lithium yield from AGB 
stars of different mass may contribute to understand the role 
(if any) of these stars in the formation of the SG in GCs. 
It is commonly believed that the polluting matter must be 
diluted with pristine matter to explain the abundance pat- 
terns, such as the Na-O anticorrelation (Prantzos & Car- 
bonnel 2005, D'Antona & Ventura 2007). If the progenitors 
of the SG stars are massive stars, for sure they do not have 
any lithium left, and the lithium in the SG must be due to 
the mixing with pristine gas. If instead the progenitors are 
massive AGB stars, they may have a non negligible lithium 
yield, that must be taken into account in the explanation of 
the SG abundances. 



2 THE PROBLEMS POSED BY 
OBSERVATIONS 

Figure [T] shows a compact summary of what we know about 
lithium abundances in the halo and in GCs in the plane 
loge(Li) versus [Fe/H] 0. The halo stars are plotted as tri- 
angles, from Melendez et al. (2009). We show both the LTE 
and non LTE abundances (these latter as grey triangles). 
The data for th ree clusters are added, at their [Fe/H] con- 
tent, taken from lCarretta et al.l ()2009d ) scale. The references 
for the clusters data are in the Figure label. Notice that the 
three open triangles of NGC 6397, at much lower e(Li) than 
the other points, refer to subgiants, in which Lithium is re- 
duced by mixing. Although the data analysis is not homoge- 
neous among the different samples, the figure shows interest- 
ing trends. The lithium spread in the field stars in the range 
of metallicities of the clusters NGC 6397 and NGC 6752 is 
very small around a plateau value loge(Li)~2.2 for the non 
LTE abundances, and slightly larger for the LTE determi- 
nations. In fact the triangles at loge( Li)< 2 are lower mas s 
stars for which depletion is expected (|Melendez et al.ll2009l ). 
The WMAP - big ba ng nucleosynthesis " standard" abun- 
dance, loge(Li)=2.72 (jCvburt et all |2009| ) is much larger 
than the plateau abundance. The li thium spread in th e clus- 
ters appears a bit larger, although iLind et al] (|2009l ) point 
out that in NGC 6397 it is consistent with the observa- 
tional error. We should expect a larger lithium spread among 
GC stars if there are SG stars, even if the pollutors' gas 
(AGB or ma ssive stars envelopes) has been diluted wit h 
pristine gas (jPecressin et al. I l2007bl : iPrantzos et al] |2007| ). 
The dilution is very plausible if there is a direct correla- 
tion between lithium and so dium abundances, as convinc- 
ingly shown in NGC 6752 (|Pasquini et al.ll2005l ). A sim- 
ilar correlation in NGC 6397 is based only on the high 
sodium a bundance of the three subgiants plotted as open 
triangles (iLind et al.l 20091). and is not co nvincing among 
the stars of 47 Tuc ( Bonifacio et al.l |2007| ) that may suf- 
fer from norma l depletion m echanisms due to their larger 
iro n content (jP'Orazi et al.i ,2009.) . In addition, according 
to iPasauini et al.l l|2008l ). two stars in NGC 6397 differ by 
~0.6dex in oxygen, but have "normal" loge(Li)~2.3: this 
is certainly not easily compatible with a simple dilution 
model, and may require that the pollutors are also impor- 
tant lithium producers. In fact, if the AGB pollutors produce 
enough lithium, a dilution model must take it into account. 
Notice that the dilution model is not so straightforward as 
we may think a priori: it will include a fraction a of matter 
with pristine Li, plus a fraction (1-a) having the Li of the 
ejecta (so, either the abundance of the AGB ejecta in the 
AGB mass range involved in the SG formation, or zero Li 
for the FRMS model). The dilution required to explain a 
given range in observed Li is different if we assign to the 
pristine Li the value loge(Li)=2.72 (see above), or the atmo- 
spheric pop. II value (~2.36), or some intermediate value. In 
addition, if we are assuming that the uniform surface abun- 
dance of Li in pop. II is due to a depletion mechanism, also 
the abundance resulting from the dilution model must be 
decreased to take into account a similar depletion factor. 
A different interesting problem is posed by the GCs 



^ the standard notations, i.e. loge(Li)=log(Li/H)-|-12 and 
[Fe/H]=log(Fo/H)-log(Fe/H)0, were adopted 
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in which a "blue" main sequence (MS) has b een revealed 
from precise HST photometry, namely uj Cen (|Bedin et alJ 
I2OO4I ) and NGC 2808. The blue MS can only be inte rpreted 
as a very high helium ( mass fraction Y~0.38) MS (|Norrisl 
I2OO4I : iPiotto et al.ll2005l ). Actually, in NGC 28 08 three MS 
well separated each other in color are present l|Piotto et al.l 
I2OO7I ). corresponding to three main helium content values, 
and in agreement with the predictions made from the distri- 
bution of st ars in the very extended and multimodal horizon - 
tal branch (|D'Antona fc Caloi|[20o3 : lD'Antona et al.ll2005l ). 
As the blue MSs are well detached in color, all their stars 
must have helium content in a very strict range of values 
— practically a unique value of helium. This poses a prob- 
lem for the FRMS model, in which the stars are formed in 
the discs ejected by the progenitors. In this model, the SG 
helium abundance will reflect the abundance in the disc, 
that will vary from one massive star to another, or even 
will be linked to the evoluti onary stage at which th e stel- 
lar envelope matter is lost (jPecressin et al.l l2007lj ) . Con- 
sequently, stars born out of such discs will show a spread 
in helium abundances, that will result in a broadened MS 
rather than in a well detached blue MS fR cnzini 20081 ). On 
the contrary, IPumo et al.l (|2008| ) noticed that the helium 
abundances of super-AGB stars enve lopes are within the 
small range 0.36<Y<0.38 (|Siessll2007l ), and iD'Ercole et all 
l|2008l l have shown that a full chemo-hydrodynamical model 
of the cluster can provide a reasonable interpretation of the 
three MSs of NGC 2808, provided that the blue MS is formed 
directly by matter ejected from the super-AGB range, undi- 
luted with pristine gas. In the future, spectroscopic observa- 
tions of the blue MS in u Cen and NGC 2808 will provide a 
falsification of this hypothesis, e.g. by means of the oxygen 
and sodium abundance revealed. In particular lithium can 
be an important test too, as it could provide an independent 
calibration of the mass loss rate in the super-AGB phase. 
Already some observations of lithium at the turnofi^ of uj Cen 
are available (Bonifacio et al. 2009, in preparation), but it is 
not clear whether they contain stars belonging to the blue 
MS. Our computations thus have two aims: 1) to compute 
the average lithium abundance in the envelopes of different 
evolving AGB masses, and show its dependence on the mass 
loss formulation; 2) to provide the first computation of the 
super-AGB stars lithium yields, when using the "standard" 
mass loss rate formulation adopted for the massive AGB 
stars. The results presented here can be applied to model 
the lithium abundance in the SG of different GCs. 



Table 1. Lithium average abundances in the ejecta for Z=10 
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conductive opacities are taken from Poteckhin (2006, see 
the web page www.ioffe.rssi.ru/astro/conduct/), and are 
harmonically added to the radiative opacities. Details on 
the rest of the physical inputs (equation of state, convec- 
ti on, overshooting, nuc l ear re action rates) can be found 
in lVentura fc D'Antonal (|2009l ). Time-dependent mixing of 
chemicals within convective zones has been treated as a dif- 
fusive process, following the approach by Cloutman fc EoU 
(1976), solving for each chemical species the diffusive-like 
equation: 



dJCi / dXi > 



drrir 



(1) 



where D is the diffusion coefficient, for which, given the con- 
vective velocity v and the scale of mixing I, a local approxi- 
mation {D ~ !«/) is adopted. The borders of the convective 
regions are fixed according to the Schwarzschild criterion. In 
our "standard" models, mass loss is described according to 
the Blocker (1995) formulation, that extends Reimers' recipe 
to describe the steep increase of mass loss with luminosity 
as the stars climb the AGB on the HR diagram. The full 
expression, for Mira periods exceeding lOOd, is 



M = 4.83 X 10~^^r?flA/ 



(2) 



3 THE MODELS 

In this work we show the lithium evolution obtained 
in our rece nt AGB models, dis c ussed in the series of 
papers by IVentura" fc D'Antonal (l2005allbh. and whose 
yields are summarized in Table 2 of IVentura fc D'Antonal 
l|2009l ). We deal with a metallicity of Z=10"^ in aU 
the models described here. The numerical structure of 
the evolutiona r y cod e ATON is described in details in 



Ventura et al.l (Il998l) . We adopt the latest opacities by 
Ferguson et al.l l|2005l ) at temperatures lower than 10000 
K and the OPAL opacitie s in the version documented 
by llglesias fc Rogerj (|l996l ) . The mixture adopted is a- 
enhanced, with [a/Fe] = 0.4 (|Grevesse fc Sauvallfiggsj ). The 



where r^R is the free parameter entering the Reimers' (1977) 
prescription. We use r]R — 0.02, according to the calibration 
based on the luminosity function of lithium rich stars in the 
Magellanic Clouds given in Ventura et al.(2000). In this work 
we also consider an extreme value of rjR = 0.1, and models 
adopting the Vassiliadis fc Wood (1993) mass loss rate. 

For core masses above ^l.OSM© — for these models, 
above initial mass 6.3 Mq — carbo n ignites in the C-O core, 
in conditions of semi-degeneracy (iRi tossa et al.llT999l : ISiesa 
2007) and the star acquires a degenerate O-Ne core that 
can evolve through thermal pulses like an AGB star. We fol- 
low the lithium phase in the super-AGB models of masses 
M^6.5Mq. The details of these models are presented else- 
where (Ventura 2009, in preparation). 
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Figure 2. Prom bottom to top panel we plot the luminosity, 
surface lithium, HBB temperature and ^He surface content along 
the AGB evolution of the 5 Mq with standard ?7fl=0.02 mass loss 
rate. The two vertical lines delimit the time of maximum lithium 
production, when the surface lithium abundance roughly exceeds 
loge(Li)=3. 



Figure 3. Temperature at the bottom of the convective layer 
(bottom) and lithium surface abundance (top) as a function of 
the time for the masses 6, 5 and 4Mq, from left to right. Time 
is computed from the beginning of the AGB phase, when the 
H-shoU burning is reignited. 



4 RESULTS 

We show in Figure[2]the time evolution of luminosity, surface 
lithium abundance, temperature at the bottom of the con- 
vective envelope (Tbco) and '^He surface abundance along 
the evolution of the standard (?7_r=0.02) model of SAf©. 
The time variable is set to zero at the end of the core he- 
lium burning phase. We see that the HBB phase begins 
when Tbcc~ 40 x lO^K, but the lithium-rich phase, when 
loge(Li) exceeds 3, occurs when Tbce ^ 60 x lO^K, and it 
ends ~ 2 X 10*yr after, due to the exhaustion of the ^He in 
the envelope. 

Figure [3] displays the lithium production and Tbcc in 
the envelopes of stars of 6, 5 and 4Mq. We adopt a lin- 
ear scale for the mass fraction of lithium (Xli) that shows 
more clearly the timescale of production and destruction. 
The larger masses achieve a larger peak of abundance, but 
their production timescale is shortefl This is due to the 
larger Tbcc, that causes a faster consumption of the ''He 
buffer. The lithium yield will depend on the lithium abun- 
dance reached in the envelope, on the duration and on the 
mass lost during the phase of lithium production. In Ta- 
ble 1 we show the average abundance in the whole ejected 
envelopes (the yield is obtained by multiplying the abun- 
dance to the total envelope mass, equal to the difference 
between the initial mass and the core mass listed in the 
third column). For 4, 5 and 6Mq we also show the result 
for the three different mass loss rate formulations. In the ex- 
treme case of mass loss {riu=0.1) shown in the Table, we see 
that average abundances 0.6-0.7dex larger than the stan- 
dard one s are reached, while the a bundances obtained by 
adopting IVassihadis fc WoodI l|l993l) mass loss formulation 
are negligible. The comparison of the lithium production 

^ Notice that also the total duration of the thermal pulse phase 
is shorter, for larger masses, due to the faster mass loss. 
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Figure 4. The surface lithium abundance as a function of the 
stellar mass, during the evolution with mass loss starting at 5 Mq . 
The standard case 77^=0.02 is compared to the larger mass loss 
rate case. 

during the 5Mq standard and rju— 0.1 models is shown in 
FigureUJ the abundances for these models are shown as open 
circles in Figure [T] The standard models (77_r=0.02) show a 
non monotonic behaviour: the yield from the 4MQis larger 
than at 5 Mq , in spite of the lower peak abundance reached 
in the envelope, due to the longer duration of the lithium 
rich phase. However, the 6 Mq yield is also larger than the 
5 Mq yield, thanks both to the larger abundances and to the 
larger luminosity of the star, that favours a larger mass loss 
rate. The evolution of the 6 Mq is also characterized by the 
fact that the HBB and the lithium production begin before 
the thermal pulse phase. Table 1 also includes the standard 
results with mass spacings of 0.5 M©. The behaviour of the 
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6 6.5 7 7.5 



Figure 5. Temperature at the bottom of the convective layer 
(top) and Hthium surface abundance (bottom) as a function of 
the total mass of the mass losing stars of 6, 6.5, 7, 7.5 and 8 Mq. 



that the lithium yield becomes very large in the super-AGB 
evolution: it may become much larger than the population I 
abundance for progenitor masses of ~8 Mq. This star has an 
0-Ne core of 1.36 Mq, so it is close to the mass limit of the 
super-AGB range. Of course, this very high lithium produc- 
tion is linked to our own modelling of the physics of massive 
AGB stars and super-AGB stars, and in particular to the 
choice of the mass loss law. More generally, these models pre- 
dict that the stars belonging to the blue main sequence of uu 
Gen and NGC 2808 can be richer in lithium than the other 
SG stars in GCs, if they were born directly from super-AGB 
ejecta, as we propose to explain their high and uniform he- 
lium abundance. If this prediction will be demonstrated to 
be incorrect, the mass loss rates from super-AGB stars (and 
massive AGB stars as well) must be revised downwards. As 
a consequence, the lithium production from massive AGB 
stars will not be significant, and both dilution models, either 
with AGB stars or with FRMS matter, will not provide dif- 
ferent informations for lithium. We then urge that a strong 
observational effort is made to observe spectroscopically the 
blue MS stars. 



lithium average abundance is not perfectly monotonic, but 
there is a general trend (decreasing abundances for increas- 
ing initial mass) that is inverted as soon as we deal with the 
super-AGB models. The physical conditions in super-AGB 
stars are in fact more extreme: their luminosity is larger, 
so that we have to deal with larger mass loss rates, and 
they achieve a larger Tbcc, so that more extreme HBB nu- 
cleosynthesis is expected. HBB begins very early during the 
star evolution. Increasing the mass, the entire lithium pro- 
duction and destruction phase occurs before the onset of 
thermal pulses, as we show in Figure [S] Given the chosen 
mass loss law, the 8 Mq evolution actually may not even 
reach the thermal pulse phase. The yield of the 8 Mq be- 
comes enormous. Of course, a different mass loss law may 
dramatically decrease the lithium average abundances. 



5 CONCLUSIONS. LITHIUM IN THE SECOND 
GENERATION OF GALACTIC GLOBULAR 
CLUSTERS 

Our results show that the lithium yield from massive AGB 
stars is dramatically dependent on the adopted mass loss 
formulation. If we take our "standard" results in Table 1 
at face value, ignoring the big question mark on mass loss, 
the yields can be used to predict the lithium expected in 
the SG, if the SG is a result of star formation from AGB 
ejecta diluted with pristine gas. The abundances will de- 
pend mainly on the mass range of the AGB progenitors: if 
the ejecta of masses in the range 4.5-6 Mq are involved, 
their abundance is loge(Li)~2, 0.7dex smaller than the Big 
Bang abundance. In order to explain the abundances ob- 
served in NGC 6397 or in NGC 6752, a dilution model in- 
cluding the ejecta of these AGB stars will require only a 
slightly smaller percentage of pristine matter, with respect 
to a model including the lithium free FRMS, and we can 
not discriminate between the two models. The situation is 
different for the cases in which super-AGB ejecta may be 
involved. In fact the main result of these computations is 
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